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^ ■ Abstract 

O 

We investigate the one-dimensional(lD) d-p model, simulating a Cu-0 linear 

^ I chain with strong Coulomb repulsion, by using the numerical diagonalization 

OO _ 

'sj" I method. Using the Luttinger liquid theory, we obtained phase diagrams of 

O 

o 

Q>^ I action at d-sites and the nearest-neighbor interaction between p- and d-sites 



the ground state on Ud — Upd plane, where Ud and Upd represent on-site inter- 



"i , perconducting phase (SC(I)) is restricted to attractive interaction Upd < 0. 

c ; 

Q , On the other hand, in the strong coupling region, we found a novel supercon 

ducting phase(SC(II)) for repulsive interaction Upd > \Ud\ and a insulating 
rS ' state with a charge gap for Ud > U^ and Upd > U^^ with critical values U^ 

and C/p^ at half-filling. Away from half-filling, another superconducting phase 



(SC(III)) appears for Ud >> Upd > 0; which has been found for C/^ — > cx) in 
the previous paper [Physica C205 (1993) 170]. An analysis of the spin gap 
suggests that the SC(I) and SC(II) include the Luther-Emery region (with 
spin-gap) in part, while the SC(III) belongs to Tomonaga-Luttinger region 
(without spin-gap) in whole. 



I. INTRODUCTION 

Since the discovery of copper oxides superconductors, strongly correlated electron sys- 
tems have been extensively studied. In particular, there is much theoretical interest in the 
electronic state of the d-p model because of the possible relevance of high- temperature 
superconductivity [1-12]. In the previous work, the present authors have studied the one- 
dimensional d-p model, simulating a Cu-0 linear chain with infinite intra-site interaction 
f/^at the Cu-(d-)sites and the nearest neighbor interaction Upd between the 0-(p-)sites and 
d-sites, is solved exactly by using the numerical diagonalization method |jl3|. By assum- 
ing the Luttinger liquid relations, the superconducting correlation is found to be dominant 
compared with the CDW and SDW correlations in the proximity of the phase boundary 
towards the phase separation. Recently, A.Sudb0 et al. showed that the d-p chain with 
large Upd exhibits flux quantization with charge 2e and slow algebraic decay of the singlet 



superconducting correlation function on oxygen sites |jT^. These works demonstrate that the 
numerical diagonalization studies of finite sizes systems combined the Luttinger liquid the- 
ory have supplied us with unambiguous and important information about the complicated 



electronic systems such as the d-p model [13-Q 



In this paper, we study the one- dimensional d-p model by the numerical diagonalization 
method. To clarify the role of the inter- and intra-site interaction for the electronic structure 
of d-p model, we turn our attention to Ud and Upd- In the weak coupling regime, we will 
take advantage of the weak coupling theory (so-called g-ology) to analyze numerical results. 
For Upd -^ cxD, we find some exact results which agree with numerical results. This paper is 
organized as follows. In the next section, we define the model Hamiltonian. The Luttinger 
liquid relation is also discussed. In Sec. 3, we present our numerical results of phase diagram 
on the Ud — Upd plane by using the Luttinger Liquid relation. We also analyze a spin gap. 
The conclusion of this work is presented in Sec. 4. 



II. MODEL AND LUTTINGER LIQUID RELATION 

We consider the following model Hamiltonian for the Cu-0 chain: 

H = -t J2 (Ptdj'^ + ^-c-) + ed ^ d]^dj„ + tpj^ ptPia 

<ij>,cr j,cr i,a 

+ UdY^ ndj^ridji + Upd Yl ^pi^ndja' , (1) 

j <ij>,<ja' 

where dj^ and p]^ stand for creation operators of a electron with spin a in the d-orbital at 
site j and of a electron with spin a in the p-orbital at site i, respectively. Udja = dj^dj^ and 
^pia = PiaPicr- t stauds for the transfer energy between the nearest neighbor d- and p-sites, 
which will be set to be unity (t=l) hereafter in the present study. The atomic energy levels 
of p- and d-orbitals are given by ep and e^, respectively. The charge-transfer energy A is 
defined as A = e^ — e^. 

The fillingness n is defined hj n = N^/Nu, where A^„ is the total number of unit cells 
(each unit cell contains a d- and a p- orbital), N^. is the total electron number and the Fermi 



wave number kp is given as kp = f n. In the non-interacting limit with Ud = Upd = 0, the 
hybridized bands are given as 



E^it) ^ ""'-"f "^'l (2) 

where tk = 2tcos{k/2) and k = 27ri/Nu {i = 0, ±1, ±2, ...)• Here, E+{k) and E-{k) stand 
p- and d- hke bands respectively. 

To achieve systematic calculation, we use the periodic boundary condition for A^e = 
Am + 2 and antiperiodic boundary condition for Ne = 4m with an integer m. This choice of 
the boundary condition removes accidental degeneracies so that the ground state is always 
a singlet with zero momentum. We numerically diagonalize the Hamiltonian with up to 12 
sites (6unit cells) using the standard Lanczos algorithm. We have calculated the ground 
state energies per unit cell for different sizes of systems. The relative difference is of order 
of 10~^ ~ 10~^ for typical values of parameters involved in the model. It indicates that size 
dependence in the ground state energies is negligible. 

The chemical potential fi{Ne, N^) is defined by 

,iN,_,N^) = Mf,±MA^Mf,^MA , (3) 

where EQ{Ne, Nu) is the ground state energy of a system with A^„ unit cells and N^ electrons. 
When the charge gap vanishes in the thermodynamic limit, the uniform charge susceptibility 
Xc is obtained from 

Xc{Ne, Nu) = ... ^. ... ^TTT -^rrr ■ (4) 

H{Ne + 1, Nu) - H{Ne - 1, Nu) 

We will discuss the correlation functions in connection with the Luttinger liquid the- 



ory p!8|-p0|. Some relations obtained in the Luttinger liquid theory have been established 
as universal relations in various one-dimensional models. Some of one-dimensional models 
can be solved rigorously by the combined use of the Bethe ansatz method with the nu- 



merical methods and/or the conformal filed theory [pi|-pB|]. In these theories, the critical 
exponents describing the power-law decay of various types of correlation functions have 
been determined. The bosonization theory on the Tomonaga-Luttinger model p4| , p5| and 



the Luther-Emery model |^6[ and the week coupling renormalization group theory (known 



as g-ology) also provides us with some rigorous information about the critical exponents 



27| , |28[| . In the Luttinger liquid theory, the critical exponents of correlation functions are 
determined by a single parameter Kp for isotropic models in spin space, or more explicitly, 
some correlation functions in momentum space have singularities as follows: 

1) Tomonaga-Luttinger (T-L) regime (without spin-gap) 

CcDw{k) and CsDw{k) r^ \k - 2kF\^'' for k r^ 2kF , (5. a) 

\k\^ for fc ~ , 

Cssik) ~ <! ' ' , ^ ^ (5.6) 

\k-2kF\^ " for kr^2kF , 



Crsik) r^ \k\^p for A; ~ , (5.6) 

2) Luther-Emery (L-E) regime (with spin-gap) 

CcDwik) r^ \k - 2kF\^''-^ for A; ~ 2A;j. , (6.a) 

Css{k) ^ \k\^'^ for fc~0 , (6.6) 

Cs'Dvy(A;) and Crsik) no singularity , (6.c) 

where C^Dvy; C'sDvy, C'ss and Cts stand for Charge Density Wave (CDW), Spin Density 
Wave (SDW), Singlet Superconducting and Triplet Superconducting correlation functions 
respectively. In the Luther-Emery regime, the spin excitation spectrum has a gap, while in 
the Tomonaga-Luttinger regime, the spin is gapless. In the T-L regime, spin excitation, the 
SS and TS correlation have the same critical exponent apart from a logarithmic correction. 
On the other hand, in the L-E regime, the real space correlation functions of SDW and 
TS decrease exponentially with distance r. The parameter Kp is related to the charge 



susceptibility Xc and the charge velocity Vc by the relations |18,20,22 



71" 

^P = 2^cXc , (7) 

v,= -^(E^~Eo) , (8) 

where Ei — Eq is the lowest charge excitation energy. These relations tell us that the enhance- 
ment of Xc leads to increase of absolute value of the critical exponent in the superconducting 
correlations and decrease of those in the CDW and the SDW correlations. 

For the one- dimensional d-p model in the weak coupling regime, T. Matsunami and 



M. Kimura [^ have calculated critical exponents in some correlation functions by the 
renormalization group analysis. They showed that the d-p model is mapped onto ^f-model 
with spin independent couplings : 

9i = Ud\atJ^ + AUpd cos {kF)\akJ'^\atJ'^ (9.a) 



where la^^P = |(1 ± A/wA^ + 4t|). The T-L regime corresponds to the case with gi > 0, 
while the L-E regime corresponds to the case with gi < 0. In the most divergent approxima- 
tion, the critical exponent Kp is given hy Kp = 1 + {gi — 2g2)/{2nvp). The superconducting 
phase appears in the region of gi — 2g2 > 0. The transformed phase diagram on the Ud-Upd 
plane is shown in Fig.l. 



III. NUMERICAL RESULTS 

In the non-interacting case, the critical exponent Kp is unity at any A and any filling for 
infinite systems. Numerical results of Kp indicate that Kp is equal to 0.92 for the system of 



6- unit cells with 8-electrons and 0.90 for 4- unit cells with 4-electrons at A = 2. It suggests 
that the size dependence of Kp are not so large even the 4-units system. 

For finite coupling case, the numerical values of Kp are 0.86 (0.91) ,0.80 (0.81), and 
0.77 (0.72) for Ud =0.5, 1.0, and 1.5 respectively in the system of 6-unite cells with 8- 
electrons at A = 2 and Upd = 0, where the values in parentheses are obtained by the 
g-ology. The consistency of both results is good in the weak coupling regime. We expect 
that our numerical results are sufficiently reliable for the strong coupling regime. 



A. Phase diagram at half-filling 

First we investigate the electronic state at half-filling (n=l). In Fig. 2, we show the 
contour map for Kp on the Ud — Upd plane at A = 2 for the system of 6-unit cells with 6- 
electrons, where the contour lines are drawn by using the spline interpolation. The numerical 
result agrees with the result of g-o\ogj in the weak coupling regime. It indicates that the 
superconducting phase (SC(I)) appears in the lower half-plane ( Upd<0). The attractive 
interaction Upd accounts for this superconducting phase. Moreover, in the strong repulsive 
regime, a novel superconducting phase (SC(II)) appears for the large repulsive interaction 
Upd in the proximity of the phase separation (PS (II)). These results seem to be similar to 
that of the one band model ( U- V model) obtained in the previous work []l6i|T7 . 



We also evaluate the charge excitation gap Ac by the discontinuity in the chemical 
potential at half-filling {n = 1) and show the region of the insulating state with finite Ac in 
Fig. 2. For large Ud and Upd, we find that the charge gap exists on a large area. 

In the case oi Ud = oo, we have shown that the charge gap Ac is roughly proportional 
to A for large A |13[ . By using the least square method, we estimate Ac in more detail as a 



function of A for Ud = oo and Upd=0 in Fig.3. It indicates that Ac is roughly proportional 
to A — 2t^/A for large values of A |^. The effect of Upd is to increase Ac by roughly 2Upd, 



which corresponds to the mean field contribution of Upd to Ac. Accounting this effect, we 
estimate the critical value Up^ of the metal-insulator transition at A = 2. In this case, we 
get Ac ~ 1.3 for Upd = and f/^^ = —Ac/2 ~ —0.7. This result consists with the numerical 
result in Fig. 2. 

On the other hand, in the limit Upd —>■ C)0, the problem becomes very simple and we can 
get some exact results. The exact wave function of the insulating state is given by 

\ins >= I ,1,0,1,0,1,0,1,0,1,0, >, 

where each underlined number denotes the number of electrons at a Cu-site and another 
does the number of electrons at a 0-site. The energy of this insulating state is Eins = 0. We 
also consider a metallic state defined by the first charge- excited state from the insulating 
state, whose element is presented by 



\metal >= |...., 0, 1, 0, 2, 0,h, 0, 1, 0, > 



where "2" denotes a doubly occupied site ( we call this electron pair dimer ) and "/i" 
denotes a hole which can move freely in the d-like band. The energy of the metallic state 



is E^etai = Ud + E-{k = 0), where E-{k) is defined by eq.(2) E-{k = 0) = -1.236.. in the 
case of A = 2. When Emetai > E^s, the insulating state is the ground state. Then, the 
metal-insulator transition occurs at U^ = —E~{k = 0) in the limit Upd — * oo. Our numerical 
result agrees with the exact result mentioned above. 

When Ud is smaller than f/J, many charge excitation occur; which make dimers and holes. 
To gain the kinetic energy of holes, dimers are required to aggregate each other. Therefore, 
the phase separation (PS (II)) occurs such as 

|PS(II) >= |...,2,0,2,0,2,0,2,0,l,0,h,0,l,0,h,0,l,0,... > 

as seen in Fig. 2. When Ud becomes large negative value, all of the electron make dimers and 
the phase separation disappears. 



B. Phase diagram away from half-filling 

Next we investigate the electronic state away from half-filling. In Fig.4, we show the 
phase diagram on the Ud — Upd plane for n = 4/3( 8-electrons per 6-unit cells). 

In the attractive coupling regime with Ud < or Upd < 0, the phase diagram is similar 
to that of half-filling. In the repulsive coupling regime with Ud > and Upd > 0,the phase 
separation region (PS (III)) appears instead of the insulator phase seen in half-filling. Fur- 
thermore, the new superconducting phase (SC(III)) is also found in the proximate region 
towards the phase separation. The mechanism of the phase separation PS (II) can be un- 
derstood as well as that of the half-filling case. However, for PS (III), the mechanism of the 
phase separation changes completely, since double occupancy of electrons on each Cu-site 
is prohibited and doped electrons (they make dimers) must sit at 0-sites at the cost of the 
atomic energy difference A. In the limit of Ud and Upd — !> oo, a schematic wave function of 
PS (III) is given as 

|PS(III) >= |...., 0, 2, 0, 2, 0, 2, 0, 0, 1, 0, 1, 0, 1, 0, .... > . 

Again, dimers are aggregated each other to gain the kinetic energy. Of cause, for Ud < 2A, 
dimers can sit also on Cu-sites and the above type of the phase separation vanishes. It agree 
with our numerical result in Fig.4. 



C. Analysis of Spin gap 

Finally, we consider the energy gap in the spin excitation spectrum. For weak coupling 
regime, the boundary line between T-L regime and L- E regime obtained by the g-ology on 
the Ud — Upd plane is gi = 0, i.e., Ud = —iUpdCos (/cF)|afcpP/|afcpP (See Fig.l). For a fixed 
value of Upd, the spin gap vanishes at a critical value t/J. In the weak coupling regime, C/J 
increases with Upd for n < l{kp < vr/2). On the other hand, in the strong coupling limit 
Upd — ^ oo,the spin gap vanishes at the critical value f/J = 2E~{k = 0) = —2.472 exactly 
3^ . To presume the boundary line for intermediate coupling regime, we will estimate the 
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spin gap numerically. The spin excitation energy for a finite size system is obtained from the 
energy difference between the lowest triplet state and the singlet ground state. We assume 
that the size dependence of the spin gap as 

As{N^Y = AsiooY + C/Nl (10) 

where As(A^u) is the spin gap of A^^j-site system and C is a constant. Sano and Takano ||^ 
have shown that this finite size scaling is successful in the estimation of the spin gap for 
t — J — J' model. Hereafter, we denote the As(oo) as As for simplicity. 

We calculated the As as a function of Ud under fixed values of Upd using the eq.(8) with 
the Nu =3 and 6 systems for n = 4/3. Fig. 5 shows that a large spin gap opens in the 
region Ud< — 2. It decreases with increasing Ud and seems to close to zero at f/^ ~ for 
Upd = 0. We get f/j> — 0.5 for Upd =1,2 and 3 and f/j ~ —1 for Upd =5. Because of the 
finite size effect, it is difficult to determine the ?7J accurately. The result of the limiting case 
Upd -^ C)0 (f/j = —2.472...), the g-o\gy and the numerical data for Upd < 5 suggest that the 
superconducting phases SC(I) and SC(II) belong to L-E regime in part, while the another 



superconducting phase SC(III) belong to L-E regime in whole ^3 



IV. SUMMARY AND DISCUSSION 

In summary, we have numerically diagonalized the one-dimensional d-p model with finite 
sizes. Paying special attention to the role of Upd in the strong coupling regime, we have 
calculated the critical exponent of correlation functions, charge gap and spin gap in the 
systems. On the Ud — Upd plane, the obtained phase diagram consists with the result of 
g-ology in the weak coupling regime. On the other hand, in the strong coupling regime, 
more complicated phases are obtained. Using the Luttinger liquid relations, we found the 
three superconducting phases (SC(I),SC(II) and SC(III)) in the proximity of the phase 
separation region ( PS(I),PS(II) and PS(III) ) respectively. The relation (8) tells us that 
the enhancement of Xc leads to increase of absolute value of the critical exponent in the 
superconducting correlations and decrease of those in the CDW and the SDW correlations. 
Thus, unstable phase is accompanied by the superconducting phase. 

We speculate that the three types of mechanism of phase separation produce the three 
types of superconducting phases. The origin of unstable phase PS (I) and the supercon- 
ducting phase SC(I) are explained by the attractive interaction Ud and/or Upd- The large 
repulsion Upd produces the unstable phase PS(II) and the superconducting phase SC(II) at 
half-filling (similar phase has been obtained in the ID U — V model) [26|,|l7|. Away from 
half-filling, the third unstable phase PS(III) and the superconducting phase SC(III) appear 
in the region ^7^ >> and Upd > . Our result suggests that Upd enhances the charge 
fiuctuation which promotes the superconducting correlation. It is noted that the parameter 
region of SC(III) is corresponding to the charge transfer insulator region at half-filling. We 
believe that this result have the relevance to the high-temperature superconductivity which 
is realized by doping to the charge transfer insulator. 

The spin gap analysis suggests that the SC(I) and SC(II) phases partially include the 
Luther- Emery region (with spin-gap). At the same time, the whole SC(III) phase belongs 



to Tomonaga-Luttinger liquid (without spin-gap). Generally, in realistic parameter region 
Ud » Upd > 0, the system belongs to the T-L regime. It is also checked as follows, we 
investigate the finite size correction to the ground state energy Eq. The conformal filed 
theory predicted the size dependence of Eq as 

Eo/Nu = e-- — + o{N^ ), (11) 

u 

where Vg is the spin velocity and c is the central charge. If the system is a Luttinger liquid, 
we should have c=l. This was shown explicitly for ID Hubbard model as well as for the ID- 
supersymmetric t — J model. For ID d—p model, we calculate the central charge numerically 
by fitting the ground state energy to the formula of the above equation by taking n =4/3 
with the system size A^„=3, 6. The charge and spin velocity, v^. and Vg are estimated for the 
largest 6-unit cluster. It turned out that c is almost constant (0.9 < c < 1.1) in the region 
Upd = 0.0 ~ 11.5 for A = 4 and Ud = 00. It suggests that the ID d — p model behaves as 
the Tomonaga-Luttinger liquid for Ud » 0. 
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Figure Captions 

Fig.l The phase diagram on the Ud — Upd plane obtained by the g-ology, where gi is given 
by eq.(lO.a) and gi — 2g2 = —{Udlatpl"^ + '^Upd\a'i:^\^\a'U'^{2 — cos(/cf))). Note that, 
in the repulsive regime Ud > and Upd > 0, {gi — 2g2) is always negative and then the 
superconducting phase does not appear within the f^-olgy. 

Fig. 2 Contour map for the Kp on the Ud — Upd plane at A = 2. We used the 6-unit 6- 
electron system (^=1) and calculated the values of Kp at points of Ud = —5, —4, ....,10 
for Upd = —2,-1, ...,10. The contour lines are plotted by using the spline interpolation. 

Fig. 3 The charge gap Ac as a function of A for n = 1. A dotted line represents A — 2t^/A, 
which is obtained the second order perturbation method (See Ref. |3( 



Fig. 4 Contour map for the Kp on the Ud — Upd plane for n=4/3 at A = 2. We used the 
6-unit 8-electron system. 

Fig. 5 The spin gap As as a function of Ud for [7^^=0,1,2,3 and 50 at n =4/3. 
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